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Abstract: 

Properties of the hadronic final state in photoproduction events witli large trans- 
verse energy are studied at the electron-proton collider HERA. Distributions of the 
transverse energy, jets and underlying event energy are compared to pp data and 
QCD calculations. The comparisons show that the "yp events can be consistently 
described by QCD models including - in addition to the primary hard scattering 
process - interactions between the two beam remnants. The differential jet cross 
sections da/dE^^ and da/drj^'^^ are measured. 
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1 Introduction 



The interaction of almost real photons and protons is the dominant process at the electron- 
proton collider HERA. A small fraction of these events have large transverse energy Et in 
the hadronic final state, measured with respect to the electron-proton beam axis, or show the 
formation of jets. 

The main characteristics of these events can be described by perturbative QCD calculations 
which are usually based on one hard parton-parton scattering per event. Two kinds of processes 
contribute to the high-E'T^ jet production in photon-proton collisions (Fig. |l|): (i) direct photon 
processes, where the photon couples directly to a parton in the proton, and (ii) resolved photon 
processes, in which the scattering occurs between a parton from the photon and one in the 
proton. Predictions for cross sections are obtained by the convolution of the matrix elements for 
the parton scattering with the parton distributions in the photon and the proton. Comparisons 
of measured cross sections with QCD calculations can provide important information on the 
parton scattering processes Q. 

Comparisons of HERA data and QCD calculations revealed that the observed jets are not 
well described by such calculations, even if they include phenomenological models for the QCD 
radiation effects and the fragmentation phase. The energy flow adjacent to the jets - the so 
called underlying event energy, or jet pedestal - was found to be far above the QCD expectations 
[§, ^, ||. This underlying event energy is measured outside of jets and includes energy resulting 
from radiation effects of the hard scattered partons and energy from the fragmentation process 
of these partons and the two beam remnants, i.e. the spectator partons. If interactions between 
the spectator partons occur in addition to the hard parton scattering, they also add to the energy 
level of the underlying event. If this underlying event energy is integrated into the measured jet 
energy, it alters drastically the measurement of the jet production rate. These effects have to 
be understood before drawing conclusions on the parton processes based on the comparison of 
measured jet cross sections with QCD calculations. 

An excess of the measured underlying event energy above the QCD calculations has previ- 
ously been observed in high energy pp collisions Q. The data could be described by adding 
interactions of the beam remnants to the calculations Q . Since the dynamics of hadronic final 
state production in photon-proton interactions is expected to be similar to that in hadron-hadron 
collisions, it is interesting to test such models in high energy photoproduction at HERA [^, ^]. 




Figure 1: Examples of diagrams for direct (a) and resolved photon (b) processes in electron- proton 
scattering. The resolved photon interactions can have - in addition to the parton scattering process shown 
- interactions of the spectator partons of the photon and the proton. 



The photon offers a unique probe for the study of the underlying event energy. Resolved 
photoproduction events can in principle involve interactions between spectator partons of the 
proton and the photon. They can be compared to direct photoproduction events, which by 
definition cannot have these additional interactions. Such comparisons give information on the 
underlying event, as will be shown below. 

In this paper detailed studies of photoproduction events with large transverse energy are 
presented. Three topics are considered: 

1. The Ej- distribution: differential cross sections da/dEj- and energy flow in the 7p center- 
of-mass system (cms) are analysed, as well as the energy fraction which is contained in 
jets. 

2. Properties of the underlying event: the energy density outside the jets is determined 
depending on the fraction of the photon energy which is available to the hard scattering 
process. Correlations of the transverse energy are measured. 

3. Jet shape and rate: properties of jets are analysed depending on the jet energy. Jet cross 
sections are presented and the influence of the underlying event on the cross sections is 
studied. 



2 Detector description and event selection 

A detailed description of the HI apparatus can be found elsewhere The following briefly 
describes the components of the detector relevant to this analysis, which makes use of the 
calorimeters, the luminosity system and the central tracking detector. 



The liquid Argon (LAr) calorimeter [10| extends over the polar angular range 4° < 6 < 153° 
with full azimuthal coverage, where is defined with respect to the proton beam direction 
(+z axis). The calorimeter consists of an electromagnetic section with lead absorbers, corre- 
sponding to a depth of between 20 and 30 radiation lengths, and a hadronic section with steel 
absorbers. The total depth of the LAr calorimeter varies between 4.5 and 8 hadronic interaction 
lengths. The calorimeter is highly segmented in both sections with a total of around 45000 cells. 
Test beam measurements of LAr calorimeter modules have demonstrated energy resolutions of 



(t{E)/E Sri 0.12/V^ © 0.01 with E in GeV for electrons ^ and cr{E)/E w 0.5/ VE © 0.02 for 



charged pions |12|. The hadronic energy scale and resolution have been verified from the ba- 
lance of transverse momentum between hadronic jets and the scattered electron in deep inelastic 
scattering events and are known to a precision of 5% and 10% respectively. 

The backward electromagnetic calorimeter (BEMC), with a thickness of 22.5 radiation 
lengths, covers the region between 151° < 6 < 177°. A resolution of 0.10/\/!E © 0.03 with 
E in GeV has been achieved for electrons. The BEMC energy scale is known to an accuracy of 
1.7%. 

The calorimeter is surrounded by a superconducting solenoid providing a uniform magnetic 
field of 1.15 T parallel to the beam axis in the tracking region. Charged particle tracks are 
measured in two concentric jet drift chamber modules (CJC), covering the polar angular range 
15° < 9 < 165°, and a forward tracking detector (FTD), covering the range 7° < 6 < 25°. The 
CJC is interleaved with inner and outer layers of multi-wire proportional chambers, which were 
used in the trigger to select events with charged tracks pointing to the interaction region. 



The luminosity system consists of two TlCl/TlBr crystal calorimeters having a resolution 
of a{E)/E = 0.1 /^/E with E in GeV. The electron tagger is located at z = —33 m and the 



photon detector at z = —103 m. The electron tagger accepts electrons with an energy fraction 
between 0.2 and 0.8 with respect to the beam energy and scattering angles below 6' < 5 mrad 
{9' = 7r-0). 

The events used in this analysis were taken during the 1993 running period, in which 26.7 
GeV electrons collided with 820 GeV protons in HERA. The scattered electron was measured in 
the electron detector of the luminosity system. This requirement restricts the negative squared 
four momentum transfer of the photon to below < 0.01 GeV^. Events were selected if they 
fulfilled the following criteria: 

1. The energy deposited in the electron tagger was in the range 8 < Etag < 20 GeV. This 
interval corresponds to scaled photon energies y = E^/Eg between 0.25 < y < 0.7. 

2. At least one charged particle was measured in the central tracker with transverse momen- 
tum above 0.3 GeV coming from the interaction region. 

3. The reconstructed event vertex lay within ±3(7 of the nominal interaction position along 
the beam axis. The width of the vertex distribution along the beam axis was approximately 
Gaussian with o" = 10 cm. 

Three different event samples were selected: 

a) A sample of events which fulfills the requirements 1,2,3 is called 'minimum bias sample'. 

This sample is only used for comparisons with results obtained with the following two 
samples (see sections 4.2 and 5.1). 

b) An event sample with large transverse energy in the final state (the 'high-i?T sample') was 

defined by Et > 20 GeV in the pseudo-rapidity range —0.8 < rj < 3.3. It contained 3254 
events and is used in sections 4.1, 4.2, 4.3 and 5.2. 

c) Events with at least one jet as defined below (the 'jet sample'). This sample contained 3499 

events and is used in sections 5.1, 6.1 and 6.2. 

The pseudo-rapidity in the laboratory system is calculated via r] = — In (tan (9/2)). The pseudo- 
rapidity in the photon-proton cms is calculated via r/* = r/ — 0.5ln Ep/E^, where — r/* ~ 2. 
The jet reconstruction was based on purely calorimetric measurements using a cone algorithm 
[|l3|| in a grid of cells in the azimuthal angle ip and laboratory pseudo-rapidity rj which extends 
from —2 < ry < 3. The cone radius R = •\/Arp + A(^ was chosen to be -R = 1.0. Each jet has 
transverse energy above E^* > 7 GeV and lies in the pseudo-rapidity range —1 < rj^'^^ < 2.5. 

The events of samples a) and b) were triggered by a coincidence of an electron tagger signal 
and a charged particle measured in the central proportional chambers. For all jet analyses, 
sample c), also one track from the central jet chamber trigger was required. Background from 
non-ep scattering is negligible in all samples. The integrated luminosity for events with (without) 
jet requirement corresponds to 290(117) nb^^ determined with an accuracy of ±5%. 



3 QCD generators 

For comparisons with the data three different event generators were used. They are based on 
tree-level QCD matrix elements. All calculations were done with the same parton distributions: 



for the proton structure the GRV-LO [14| leading order parton density parameterizations were 



used, and for the photon structure the GRV-LO |15] leading order parameterizations. As will 
be seen below, the generators can be grouped into those which include - in addition to the hard 
parton scattering - interactions of the beam remnants and those which do not. The models for 
the remnant interactions are based on parameterizations of results from hadron scattering, or 
(semi- ) hard scattering between spectator partons, or soft and (semi- ) hard parton scattering 
processes between spectator partons. 



3.1 PYTHIA 



The PYTHIA 5.7 event generator for photon-proton interactions ||T^ was used together with a 
generator for quasi-real photons. PYTHIA is based on leading order (LO) QCD matrix elements 
and includes initial and final state parton radiation calculated in the leading logarithm approx- 
imation. The strong coupling constant as was calculated in first order QCD using Aqcd = 200 
MeV with 4 flavours. The renormalization and factorization scales were both set to the transverse 
momentum pt of the partons emerging from the hard interaction. Since the QCD calculation 
of a hard parton scattering process is divergent for processes with small transverse momenta 
Pt, the requirement pt > p^^"- = 2 GeV was made. For the fragmentation process the LUND 
fragmentation scheme was used (JETSET 7.4 In this PYTHIA calculation the underly- 

ing event energy is generated by initial and final state parton radiation, and by fragmentation 
effects. 

Within PYTHIA, interactions additional to the primary parton-parton scattering may be 
generated, so-called 'multiple interactions' ||^ These are calculated as LO QCD processes 
between partons from the photon and proton remnants. Multiple parton scattering has been 
studied in proton-antiproton collisions before |19, 18, |2^. The PYTHIA version extends the 



concept of the hard perturbative QCD parton scattering to the low transverse momentum, or 
semi-hard interaction region. As mentioned above, in this low pt region the LO parton scattering 
cross section aparton diverges and becomes larger than the measured non-diffractive photopro- 
duction cross section cTnd- Since in the multiple interaction model each of the two incoming beam 
particles may be viewed as a beam of partons, the problem of too large parton cross sections 
can be solved by allowing for several parton scattering processes within one observable 7p event. 
To prevent a rapid rise of the jet cross section at small transverse jet energies, a unitarization 
scheme has been introduced which results in the damping of the cross section at small parton 
momenta. In this way, the calculated hadronic cross section can stay below the non-diffractive 
7P cross section. 

In the simplest version of the PYTHIA multiple interaction model used here, the transverse 
momentum cut-off of the hard interactions is lowered to p^'^"- < xhe mean number of 

(semi-) hard interactions is given by < n >= (Jparton{pT^"') / f^nd^ the fluctuations are calculated 
from a Poisson distribution. The number of additional interactions is typically of order 1 — 2. 
The parton process with the highest transverse momentum in the partonic final state can be 
given by any quark or gluon matrix element. This process includes then initial and final state 
parton radiation effects and its partons are connected to the beam remnants by strings. The 
additional parton scattering processes in the event are calculated as perturbative gluon-gluon 
scattering processes. The initial state gluon momenta of each subsequent process are related 
to the remaining energy of the beam remnants. The resulting fractional momenta are used to 
determine the parton densities of the beam remnants. 

The additional interactions contribute significantly to the transverse energy flow in the event. 
This is mainly influenced by the average number of interactions per event, the cut-off p^*'*, and 
by the structure functions used for the calculations. Using the GRV parameterizations for the 
parton distributions of the photon and the proton, an adjustment of the transverse momentum 
cut-off to = 1.2 GeV resulted in an adequate description of the energy flow next to jets 
as observed in the data (see section |5.1| ). The PYTHIA version with multiple interactions was 
used to determine corrections for apparatus inefficiencies. 



3.2 PHOJET 



The PHOJET 1.0 event generator was designed to simulate in a consistent way all components 
which contribute to the total photoproduction cross section |21|. It is based on the two- 
component Dual Parton Model [ 22 1 . The implementation of the PHOJET generator is similar 



to the Monte Carlo event generator DTUJET |23] which simulates multi-particle production in 
high energy hadron collisions. This latter generator was originally intended for the description 
of soft hadronic interactions and was then extended to hard scattering processes. In contrast to 
PYTHIA, PHOJET incorporates both, multiple soft and hard parton interactions, on the basis of 
a unitarization scheme [^. The soft hadronic processes are described by the soft, 'supercritical' 
Pomeron |25]. These processes are simulated by a two-string ansatz which allows for initial 
transverse momenta of the partons at the ends of the strings. The hard processes are calculated 
using the LO QCD matrix elements. Final state parton radiation effects are simulated using the 
JETSET 7.4 program [17|. Hard initial state parton radiation is not included in this version of 
PHOJET. The lower momentum cut-off for hard parton interactions was set to p™" = 3 GeV. 
Due to the unitarization scheme, small variations of this cut-off parameter do not have a large 
influence on the results of this generator. The model parameters which describe the soft part 
of the jp interactions have been tuned using results from proton-antiproton collisions and low 
energy photoproduction cross section measurements. For the fragmentation the LUND string 
concept is applied (JETSET 7.4 [ll^). 



3.3 HERWIG 



The HERWIG 5.8 ep generator is also based on the LO QCD calculations [26|. This program 
was designed to have as much input from perturbative QCD as possible, in order to minimize 
the free parameters. HERWIG includes a parton shower model which allows for interference 
effects between the initial and final state showers (colour coherence) [p7|]. The renormalization 
and factorization scales were set according to the transverse momentum of the scattered partons 
with a lower cut-off used around pj^™ > 2 GeV. The strong coupling constant as was calculated 
to first order using Aqcd = 180 MeV for 5 fiavours. A cluster model is used to simulate the 
hadronization effects [28|. 

HERWIG allows optionally for additional interactions of the beam remnants. These interac- 
tions are called soft underlying event and are parameterizations of experimental results on 'soft' 
hadron-hadron collisions. A tuning of the strength and frequency parameters is still in progress. 
Recently, also a model for multiple parton interactions has been developed for HERWIG |2£] 
which has not been used in this paper. 



4 Distributions of transverse energy 
4.1 Transverse energy cross section 

Cross sections as a function of the total transverse event energy Et have previously been mea- 
sured in 7p |^0| and pp scattering at different cm energies 34 ] . They fall steeply in the 



low Ej- region (soft hadronic interactions) and tend to fiatten towards large transverse energy. 
For sufficiently high Et the distributions can be equally well described by either a power law 
{Et)~" or an exponential decrease exp(— Ai^T^). 

In Fig. ^ the measured differential ep cross section da/dEx is shown. For this analysis the 
high-E'r event sample, defined in section |2|, was used. The transverse energy was integrated over 
the cms pseudo-rapidity range —2.5 < rj* < 1. This region is fully covered by the electromagnetic 
and hadronic sections of the LAr calorimeter. The measurement corresponds to scaled photon 
energies between 0.3 < y < 0.7 and negative squared momentum transfer < 0.01 GeV^. 



Corrections for apparatus inefficiencies applied to the observed transverse energy were pa- 
rameterized as a function of the pseudo-rapidity of the energy deposits. They decrease with 
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Figure 2: a) The differential transverse energy ep cross section integrated in the pseudo-rapidity 
interval —2.5 < rj* < 1. Full circles are data, the inner error bars are statistical, the outer error bars 
are the quadratic sum of statistical and systematic errors. The histograms show the calculations of QCD 
generators with interactions of the beam remnants (full=PHOJET, dashed— PYTHIA) and without them, 
(dotted^PYTHIA). 

b ) Corrected transverse energy flow per event versus the pseudo-rapidity in the jp cm,s, where rj* > 
refers to the proton direction. Full circles refer to events with total transverse energies between 25 GeV < 
Et{—2-5 < r/* < 1) < 30 GeV. Open circles refer to minimum bias data. Histograms are as in a). 



pseudo-rapidity and amount to factors between 1.2 — 0.8. The corrections depend on the Monte 
Carlo model used to determine them. This leads to a contribution of 10 — 15% to the systematic 
error, increasing with decreasing pseudo-rapidity. Corrections were apphed for the acceptance 
of the electron tagger and these contribute 5% to the systematic error. This error also includes 
the uncertainty of the luminosity measurement. The determination of the trigger efficiency has 
a 2% error. The dominant systematic error comes from the uncertainty in the knowledge of the 
LAr calorimeter energy scale. This error increases with Et and amounts to 15 — 50%. 

The differential cross section in Fig. ^a is compared to different QCD calculations: the 
PYTHIA generator without multiple interactions (dotted histogram) gives too small values of 
the cross section. The HERWIG generator gives similar low values (not shown). The PYTHIA 
calculation with multiple interactions (dashed histogram) gives too large cross section values 
at Et = 20 GeV, but is compatible with the data in the large Ej- region. The PHOJET 
generator, which also includes multiple interactions, provides a good description of the data 
(full histogram). 

The shape of the measured differential Ej- cross section can be described by a power law 
with n = 5.9 with a 2% statistical uncertainty in n. The power n is similar to that of the 
transverse energy jet cross section [Q]. Therefore the spectrum can be consistently interpreted 
by hard parton scattering processes. The data are in this Et range equally well described by 
an exponential function. An exponential fit gives the slope A = 0.21 at ^s-yp = 200 GeV with a 
5% statistical error in the fit. 



4.2 Transverse energy flow 

The peripheral scattering of two hadrons results in a transverse energy flow which can be de- 
scribed by a rapidity plateau of width 2 In (^/s/m), with final state hadrons of mass m carrying 
transverse momenta of around pt ~ 300 MeV, ^/s being the cms energy ||3^ . In Fig. ^ the aver- 
age transverse energy flow in photoproduction events is shown versus the 'jp cms pseudo-rapidity 
r]* . The distributions are corrected for detector effects. The open circles refer to minimum bias 
data, defined in section ^. These data exhibit a plateau within the pseudo-rapidity interval 
shown, as expected from peripheral scattering of a hadron-like photon and a proton. 

The full circles refer to the high-Sy event sample, defined in section |2[ with corrected total 
transverse energies, summed in the pseudo-rapidity range —2.5 < r/* < 1, between 25 GeV < 
Et < 30 GeV. The distribution is not compatible with being flat and increases towards the 
origin of the 'jp cms. 

Similar distributions have been observed in pp scattering with large transverse energy in the 



flnal state |33|. Here the Et distribution versus pseudo-rapidity is centered on the cms rapidity 
rj* = and has (e.g. at ^Spp = 630 GeV) a half width at half maximum of only one unit of 
rapidity for very large Et > 200 GeV which is to be compared with a possible total plateau 
width of Ar] ~ 13 (using the proton mass). 

The observed increase in the jp cms is, to some extent, unexpected: since there are relatively 
more partons at large momentum fraction x in the photon than in the proton, one would naively 
expect this distribution to be peaked somewhere in the photon hemisphere, i.e. at negative 
rapidities. The PYTHIA calculation without multiple interactions (dotted histogram) confirms 
this expectation. A multiple interaction mechanism, however, adds additional energy around the 
origin of the jp cms. The superposition of parton scattering processes by a PYTHIA calculation 
with multiple interactions (dashed histogram) results in a description which is compatible with 
the data. Here the maximum is closer to 77* = 0, but still shows a shift towards the photon 
hemisphere. The PHOJET calculation (full histogram), which also includes interactions of the 
beam remnants, provides a good description of the energy flow. 
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Figure 3: Observed relative rate of 0=circles, l=diamonds, 2=triangles, 3=squares jet events versus 
the observed transverse energy collected in the pseudo-rapidity range —3 < 77* < 1. The observed jet 
energies are above 7 GeV summed in a cone of size R ~ 1 for jet axes between —2.5 < rj* < 0.5. 
Symbols are data with only statistical error bars shown. The curves are calculations of a QCD generator 
including interactions of the beam remnants and a detailed simulation of the detector effects (PYTHIA). 

4.3 Jet rate 

In tree level QCD calculations two partons emerge from a hard scattering process, each frag- 
menting into a jet. The measured jet multiplicity depends, however, on the total transverse 
energy in the event and on the jet definition. Here a cone algorithm is used, with cone size 
i? = 1 and a transverse energy threshold E^*^-^. 

In Fig. ^ the relative rates of the observed jet multiplicities njet are shown as a function of the 
uncorrected transverse event energy Ej^^. These results were obtained using the high-E'r event 
sample. The visible jet energy threshold is set to Ej^l^-^ = 7 GeV and the jet axis is required to 
be in the range —2.5 < rj* < 0.5. The transverse event energy is the summed transverse energy 
deposited in the pseudo-rapidity interval —3 < rj* < 1. 

At Eq^^ = 20 GeV, 70% of the events have no jets at all, given the above definition, while at 
Ej!-^ = 35 GeV, half of the events contain two jets. The two jet configurations are expected from 
the tree level QCD picture. At £'j?* = 40 GeV more than 10% three jet events are observed. 



showing the effect of higher order processes with possible contributions from the photon remnant. 
At this there are essentiaUy no events without jets. 

The average fraction of the total visible event transverse energy which is contained in jets 
with the chosen parameters was calculated at Ej.^'^ = 20GeV to be ~ 10%. At Ej!-^ = 40GeV 
about half of the transverse event energy is contained in jets. 

The jet rates are by and large described by calculations of the PYTHIA generator with mul- 
tiple interactions, which include here a detailed simulation of detector effects (Fig. |3|). PYTHIA 
without multiple interactions gives much too large jet rates: at Ej^^ = 30 GeV the contribution 
of events without jets has already vanished, while the relative two-jet rate is two times the rate 
of the data (not shown). The latter model comparison indicates that the energy depositions 
beyond those resulting from the tree level hard parton scattering process plus parton showers 
are neither correlated with the jets nor jet like, but 'soft' energy depositions spread throughout 
the event. 



Summary: Distributions of transverse energy 

The photoproduction data with large transverse energy in the final state show signatures of hard 
scattering processes: the slowly decreasing transverse energy spectrum, the pseudo-rapidity 
distribution which is not flat, and the multi-jet production all demonstrate the scattering of 
constituents of the photon and proton. Comparisons of the data with different generators - 
which are based on LO QCD matrix elements plus parton showers - show that models which 
include interactions between the beam remnants provide considerably better descriptions of the 
data than models without such a mechanism. 



5 Energy of the underlying event 

5.1 Transverse energy density outside of jets 



In QCD inspired models various processes contribute to the transverse energy flow outside 
jets (e.g. @): 

A. ) Transverse energy of the partons that participate in the hard scattering including 1) initial 

and 2) final state radiation from these partons, 

B. ) Transverse energy from the interactions between the spectator partons which is essentially 

uncorrelated with the hard scattering process, 

C. ) Transverse energy of non- interacting spectator partons which is also essentially uncorrelated 

with the hard parton scattering process. 

In all three cases fragmentation effects have to be considered in addition. Monte Carlo generators 
which include all these components describe distributions of the hadronic final state in pp 
and 7p scattering considerably better than generators including only the effects of the hard 
parton interaction (item A) and remnant fragmentation (item C). The clearest experimental 
evidence for item B would be hard spectator interactions. Their observation has been reported 
in analyses of multi-jet events in pp and pp collisions [18, ^] , but were not confirmed by another 



pp experiment [19|. 



At HERA the photon offers a unique probe to study the underlying event: in jp interactions 
with resolved photons the energy flow outside jets should contain all three components mentioned 
above. Direct 'jp interactions, however, do not involve an initial state parton from the photon 
side, hence there is no initial state QCD radiation (item A.l). There are no spectator partons 
interactions (item B) and there is no photon remnant (item C). The direct processes should, 
however, have the same final state radiation (item A. 2) as the resolved photon interactions. The 
different contributions can therefore partly be disentangled by studying the underlying event as 
a function of, e.g., the fractional momentum of the remnant from the photon side (1 — x-y) which 
is for direct photon interactions and greater than in the case of resolved photon interactions. 

The jet event sample is used (defined in section ^), requiring at least two jets. To avoid a 
miss-interpretation of the photon remnant as a hard jet the pseudo-rapidity difference between 
the two jets had to be below Ar] < 1.2. The momentum fraction of the parton from the photon 
can be estimated from 

jets _ ^ + e 

where xi^^'^ is reconstructed using the two jets with the highest transverse energy E^r^^ in the 
event and their rapidities r/-'^*. The energy of the photon E^ is determined from the energy 
measured in the electron tagger. The resolution of xt^^*^ is between 15 — 30%, getting worse as 
^jeis increases from to 1. 

In this section the measurement of the energy density outside jets as a function of xl^^'^ is 
discussed. The transverse energy is summed in the central rapidity region of the 7p collision 
— 1 < ry* < 1. Energy deposited around the two jet axes within R = -y/A^j^ + Arp < 1.3 is 
excluded from the energy summation. The transverse energy density < Et > /{Ar]Aip) is then 
defined as the energy sum < Et > per unit area in the (r/, ip) space averaged over all events in 
the sample (Fig. ^. The distribution was corrected for detector effects. It rises towards small 
values of xi^'^*'^. 

Systematic errors originate from the following sources: the observed energy was corrected for 
detector effects, the correction used depends on the pseudo-rapidity of the energy deposit. There 
is a small dependence on the Monte Carlo model used to obtain the correction which contributes 
10% to the systematic error. Migration effects in x:^^** were determined with different generator 
models and lead to an uncertainty in the energy density at the level of 14%. These errors 
are shown in Fig. ^ In addition, the uncertainty in the knowledge of the hadronic energy 
scale of the LAr calorimeter may change the overall normalization of the energy density by 5%. 
Contributions from other sources of systematic errors are negligible. 

The long-dashed line in Fig. ^ indicates the energy density found in minimum bias events 
where xi^^*^ is not measurable. The energy density for directly interacting photons {xi^^^^ = 1) 
is found to be close to the energy density found in minimum bias events. This observation is 
consistent with a recent comparison of the measured transverse energy density in the central 
rapidity region of the ■yp collision. Here the same transverse energy density was measured in 
minimum bias data and deep inelastic scattering events (direct photons) |3^. These comparisons 
show that neither the initial state radiation effects of the parton from the proton, nor the final 
state parton radiation effects can be large. The measured energy density increases towards small 
xi^^*^ by 0.6GeV/rad to 3.5 times the minimum bias energy density. It is unlikely that such a 
large effect is caused alone by initial state radiation of the parton from the photon side. In the 
QCD inspired picture described at the beginning of this section, the natural explanation for the 
enhanced energy flow is interactions between the parton spectators (item B): the probability for 
these additional interactions increases with the energy of the photon remnant. 

Also shown in Fig. Q are the results of several different Monte Carlo generators. The dotted 
histogram refers to the PYTHIA calculation without multiple interactions. This calculation is 
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Figure 4: Corrected transverse energy density in the central rapidity region of the "fp collision ]77*| < 1 
per event and unit area in pseudo-rapidity and ip space, outside of the two jets with the highest Ei^* . The 
variable xi^*^ is a measure of the momentum fraction of the parton from the photon side. Circles are data, 
the inner error bars are statistical, the outer error bars result from the quadratic sum of statistical and 
systematic errors. An overall uncertainty in the hadronic energy scale of 5% is not shown in the figure. 
The long-dashed line indicates the energy density measured in minimum bias events. The histograms 
show the calculations of different QCD generators with interactions of the beam remnants (full=PHOJET, 
short-dashed=PYTHIA), and without them (dotted=PYTHIA). 



compatible with the measured energy density at large x^^*^, but shows only a small increase 
towards xl^^^ = 0, as one may expect from increased initial state radiation effects. The use of 
another parameterization of the photon structure function, e.g. LACl [^^, does not change 
this result (not shown). The HERWIG generator gives similar results (not shown). The QCD 
simulations without multiple interactions show that the observed increase cannot be understood 
as a kinematical bias. 

Adequate descriptions of the data are provided by calculations which include interactions 
between the photon and proton spectator partons, PYTHIA with multiple interactions (short- 
dashed histogram) and PHOJET (full histogram). Using PYTHIA with multiple interactions 
together with the LACl structure function parameterization instead of the GRV-LO set increases 
the number of multiple interactions per event due to the larger gluon density at small fractional 
momenta x^. A re-adjustment of the momentum cut-off from p™*'^ > 1.2 GeV to pj"*" > 2 GeV 
results again in an adequate description of the measured energy density. 

For the first time the underlying event energy has been measured in jet events using direct 
and resolved photon probes. The large difference of the energy density in the central rapidity 
region of the collision between direct and hadron like photons and the comparisons with 
different Monte Carlo models demonstrate that an additional source of Et is probably needed 
beyond that coming from a two parton scattering process and the according parton radiation 
effects. 



5.2 Transverse energy correlations 

Energy-energy correlations are sensitive measures of how energy is distributed over the available 
phase space. In this section we examine further the transverse event energy measured in Fig. ^ 
and in jet events (Fig. ^). A study of the rapidity correlations relative to the central rapidity 
region of the collision, where the event energy is largest, should provide information on the 
underlying event and is an important test of the models which describe the average event energy 
correctly. 

In Fig. m rapidity correlations ^l are shown with respect to the 7p cms pseudo-rapidity r/* = 
for the high-£^r event sample, defined in section ^. Energy deposits within jets are included in 
this measurement. The correlation function Vt was defined as 

Here Et is the total transverse energy measured in the range —3.1 < rj* < 1.3, and the other 
terms refer to transverse energies measured in the bins of size Ar] = 0.22. The average value of 
the transverse energy in a bin < E^.r]* > was determined from all events in the sample. 

The data distribution shows short range correlations around the 7p cms reference bin at 
r]* = 0. Long range anti-correlations are observed in the photon hemisphere. These have a 
minimum around rj* = —1.8 and disappear around rj* ~ —3. 

The dotted histogram in Fig. ^ shows the calculation of the PYTHIA generator without 
multiple interactions, including a detailed simulation of the detector effects. The correlations 
have the correct shape, but the (anti-) correlation is twice as strong as observed in the data. This 
implies that too much energy in the central rapidity region of the jp collision is correlated with 
the other energy deposited in the event. The PYTHIA version including multiple interactions 
is shown in Fig. || as dashed histogram. Here event energy is added around r/* = (Fig. ^). 
This stems from the interactions of the spectator partons and is not correlated with the major 
part of the event energy, which results from the primary hard scattering process. Although the 



-3 



H1 



data 

PYTHIA mia 

PYTHIA 



-2 



-1 








11 



* 



Figure 5: Observed rapidity correlations with respect to the central rapidity region of the jp collision 
T]* = 0. Full circles are data, only statistical error bars are shown. The dashed (dotted) histogram 
represents a calculation of the QCD generator PYTHIA with (without) interactions of the beam remnants. 



model does not give a perfect description of the data, it shows that the addition of uncorrelated 
energy to the events results not only in the correct average underlying event energy (Fig. 
but also gives the correct correlation strength (Fig. ^). The same conclusions hold for an event 
sample where a jet is explicitly required. 



Summary: Energy of the underlying event 

There are large differences in the underlying event transverse energy density in 7p collisions 
dependent on whether the photons involved are direct or hadron like. In the central rapidity 
region of the 7p collision, processes where large fractions ~ 1 of the photon energy are 
carried into the hard scattering process have a similar transverse energy density as minimum 
bias events. Processes with resolved photons x-y ~ are found to have 3.5 times the transverse 
energy of minimum bias events. This increase is much larger than expected from QCD generators 
using LO QCD matrix elements plus parton showers, but can be explained by models including 
interactions between the beam remnants. In addition, comparisons of the measured energy- 
energy correlations with different models demonstrate that the additional energy of the beam 
remnant interactions decreases the correlation strength to the level observed in the data. 

Therefore, the underlying event in photoproduction events can be consistently interpreted 
as the superposition of a hard scattering process plus interactions between the beam spectators. 
The characteristic kinematic quantity is the momentum fraction of the parton from the photon 
side. This determines the energy in the hard interaction and the energy left for interactions of 
the two beam remnants. 

6 Properties of jets and jet cross sections 
6.1 Jet shape 

The distribution of transverse energy around a jet axis shows an approximately Gaussian shape. 
In perturbative QCD the width of these jets is expected to decrease with increasing jet energy. 

In Fig. 1^ the uncorrected transverse energy flow in the azimuthal direction A(/9 around 
the jet with the highest E'jf* in events with at least two jets is shown. The jets were found 
in the photon hemisphere between —2 < r/jg^ < — 1 where detector corrections are small. The 

transverse jet energies were between 9 < Eji^^ < 11 GeV summed in a cone of size R = 1. 
The energy flow was integrated in a slice of \rj* — i]jet\ — ^ around the jet axis. The two jets 
usually deviate from a back-to-back configuration in the azimuthal angle. The azimuthal angle 
between the two jets was defined such that — vr < Aip < 0. Since the jets have a total width 
of approximately two units in the azimuthal direction (Fig. ^), the region 1 < A(f < 2 is 
essentially free of the second jet and shows the underlying event energy. 

The full curve in Fig. ^ represents an empirical 3 parameter fit of the jet profile using a 
function closely related to a Gaussian: 

f{A^) = Aexp{-{^\+b)' + b') + P (3) 

Parameter A describes the amplitude of the jet at Aip = which is mainly constrained by the 
energy measured in the central bin. P reflects the underlying event energy (jet pedestal). The 
full width at half maximum above the pedestal is then F = 2((ln2 + 6^)^/^ — 6)^. Regions 
potentially affected by the second jet, Aip < —1 and Aip > 2.2, are excluded from the fit. 




Figure 6: a) Observed transverse energy around the jet axis is shown versus the azimuthal angle Aip 
for jets in the photon hemisphere —2 < rj*^^ < —1 and 9 < E^r^*' < 11 GeV. The full circles represent the 
data with their statistical error bars. The axis of the second jet is between — tt < Atp < 0. The histogram 
is a calculation of the QCD event generator PYTHIA with multiple interactions, and the curve is a fit 
using eq.(\^). 

b) The fitted full width T at half maximum above the underlying event of a jet is shown depending on 
the scaled jet energy 2Ei^^ / ^JW^, where the ep center- of-mass energy is used following Full circles 

are HI jp data, the error bars represent the quadratic sum of the statistical and systematic errors. Open 
symbols are derived from fits to jet profiles ofpp data using the cms energies y^Spp Js^l , The full line 
represents a l/Ej'^* fit to the jp data. 



In Fig. ^3 the width T of jets from 7p interactions (full circles) is shown as a function 
of the scaled transverse jet energy 2E^^ / ^s^p. Here the ep cms energy is used, following a 
QCD analysis of jet shapes |^^. The decrease of T for jet energies between 5 < Ej^^ < 
20 GeV is clearly seen. This effect is well described by the QCD generators, although it is not 
inconceivable that other models could explain this effect. The decrease can be described by a 
1/E^^ dependence. Such a dependence is expected from a QCD analysis presented in |35]. 



Within errors, no dependence of P on E^^^ was found in this restricted pseudo-rapidity range. 
The amplitude A varies as a function of E^^ and T. Systematic errors in Fig. ^b (< 10%) 
are dominated by the choice of the Aip interval used for the fit. Fits to the jet profiles in the 
rapidity projection give compatible results within error bars. Since the jet pedestal increases 
as a function of the rapidity in the rj region considered here, more fit parameters are needed to 
describe the data. This leads to larger error bars on the fitted jet width than in the case of the 
fits to the (/9-projections of the jet profiles shown in Fig. 

In the same Fig. |6|b, the photoproduction data are compared to results from fits to jet profiles 
ofpp data at different cms energies (open symbols) |4^]. Here pseudo-rapidity projections 
of the jet profiles were used, allowing P in (|3|) to depend linearly on the rapidity. The same jet 
width as in ^p data was measured in pp data at different jet energies. Within the error bars the 
jet widths found in •jp and pp are compatible with having the same dependence on the ratio of 
the jet energy and the cms energy. 



6.2 Jet cross sections 

Differential inclusive jet cross sections have previously been measured in photoproduction events 
at HERA ^. The jet cross section da/dE^* as a function of transverse energy falls steeply 
and can be described by a power law {Ei^^)~^ . Such a falling distribution is expected from QCD 
calculations. Here the matrix elements of different parton scattering processes (quark-gluon, 
gluon-gluon etc.) are summed according to the quark and gluon distributions in the photon and 
the proton. The calculations are only modestly sensitive to the parton distributions, and reflect 
essentially features of the matrix elements. The measured cross section can be described by these 
calculations and can, therefore, be interpreted as the result of a parton scattering processes. 
The observed jets are mainly located in the photon hemisphere and cover there several units 
in pseudo-rapidity. Jet cross sections da/drj^^^ as a function of rapidity are sensitive to the 
parton distributions in the photon. These parton distributions can be extracted by comparison 
of the data with QCD calculations, using an iterative procedure. However, the underlying event 



energy, explained in section p.l| , has a large influence on the measured jet cross sections, as will 
be shown below. Before conclusions may be drawn on the photon structure, the underlying event 
has to be described by the QCD calculation (see Fig. or the jet data need to be corrected 
for underlying event effects (e.g. 0|). 

In Fig. the inclusive differential jet ep cross sections dcr/dE^* and da/dr]^^^ are shown. 
The measurements were made using the jet event sample. No corrections for the underlying 
event energy were applied. Here the pseudo-rapidity rf^^ in the laboratory system (rj — rj* = 
0.5 In (Ep/E^) ~ 2) is chosen for compatibility with previous measurements. The cross sections 
are for a scaled photon energy range 0.25 < y < 0.7 and negative squared momentum transfer 

< 0.01 GeV^. The jet cross section as a function of the corrected transverse jet energy in 
Fig. ^ was measured in two pseudo-rapidity intervals —1 < rj^^^ < 1 and —1 < rj^^* < 2 (see 
also tables ^). The jet cross section as a function of pseudo-rapidity in Fig. |7|b was measured 
for three different thresholds in the jet energy, E^j.''^ > 7, 11, 15 GeV (see also tables |, |, |). 
The previous measurement of da/drj^'^^ for jet energies above E^r^^ > 7 GeV suffered from 
a defect in the acceptance correction for the small angle electron detector and is superceded by 
this new measurement. 
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Figure 7: a) Inclusive differential jet ep cross sections for jets with transverse energy above 
> 7 GeV summed in a cone of R = 1 . Circles are HI data integrated in the pseudo-rapidity 



range — 1 < rf"^* < 2, squares refer to the region — 1 < rf^* < 1 (lowered by a factor 10). The inner error 
bars represent the statistical, the outer errors are the quadratic sum of statistical and systematic errors 
which depend on Ei^* or rj^^* . The uncertainty in the overall normalization amounts to 26% and is not 
shown in the figure. The curves show the calculations of different QCD generators with interactions of 
the beam rem,nants (full=PHOJET, dashed=PYTHIA), and without them (dotted=PYTHIA). 
b ) Differential jet cross section versus pseudo-rapidity for different thresholds in the jet transverse energy: 
circles are HI data with > 7 GeV, down-pointing triangles with Ei^* > W GeV, and up-pointing 
triangles with E-!^*' > 15 GeV. The dash- dotted curve is a calculation using PYTHIA without interactions 
of the beam remnants, but another parameterization of the photon structure function (see text). Other 
curves are described in a). 



Corrections to the energy response of the calorimeter were made as functions of the transverse 
energy and pseudo-rapidity of the jets. Using generator events the corrections e result from 
the comparison of jet rates resulting from generated hadrons and jets found after a detailed 
simulation of detector effects. The correction function e{E^^) decreases from 1.1 to 0.7 as 
Ej?^ increases between 7 < Ej^^ < 30GeV. The correction e{r]^^*) decreases from 1.7 to 0.7 
as r] increases from —1 to 2. The systematic uncertainty in the determination of the jet rate 
corrections is 10%, including a small dependence on the generator model used. The hadronic 
energy resolution is known to the 10% level. This was found to give a 10% contribution to the 
systematic error. The hadronic energy scale of different calorimeter segments is known to be 
the same to 3%. The latter results in an 15% contribution to the jet cross section error. The 
efficiency of the drift chamber trigger (see section |) varies as function of Ej^^ and between 
90% and 95%, with an uncertainty of 3%. In Fig. ^ the total error bars include the quadratic 
sum of the statistical and systematic errors mentioned so far. 

Sources of systematic errors affecting the cross section normalization are as follows: 1) The 
dominant systematic error of 25% results from the uncertainty of ±5% in the LAr energy scale. 
2) The error in the luminosity measurement, which includes the uncertainty in the electron 
trigger efficiency, is 5%. These errors have been added in quadrature and give a 26% overall 
systematic error in addition to the errors shown in Fig. |^. 

The fit to a power law (E'jf*)"" of the jet cross section as a function of transverse energy in 
Fig. ^ gives n = 6.1 it 0.5 for the pseudo-rapidity interval — 1 < rf^^ < 2, where the error reflect 
the statistical and systematic uncertainties of the fit. The measured power n is compatible with 
previous results on jet production in jp interactions |2|, |3|. This result compares also well with 



jet cross sections measured in pp collisions at the same cms energy -^/s^ = 200 GeV |Q where 

the same fit function, applied to the same E^^ range, results in the power n = 5.8. Within the 
errors no differences between jp and pp are found in this kinematic range. 

In Fig. |7| the calculations of the PYTHIA generator without multiple interactions (dotted 
curves) are compared with the measured jet cross sections. This model, and also HERWIG 
(not shown), show clear deficiencies in describing the cross sections of jet production at large 
rapidities and small transverse energies of the jets. This kinematic region corresponds to events 
where the momentum fraction of the parton from the photon is small Xy ~ 0.1. 

The PYTHIA calculation with multiple interactions provides a fair description of the shape 
of the measured pseudo-rapidity cross section, but gives a too large cross section in the re- 
gion of small Ef'* > 7 GeV (dashed curves). The full curves represent the calculation of the 
PHOJET generator which gives a good description of both the shape and the measured rates. 
The results of the calculations including additional event energy - PHOJET and PYTHIA - 
differ in absolute numbers in the small E^* region despite the use of the same QCD matrix 
elements and structure functions. However, the calculations include different modeling of 1) 
the beam remnant interactions ((semi- ) hard interactions in PYTHIA versus soft and (semi- ) 
hard interactions in PHOJET), 2) the transition to the non-perturbative soft scattering region 
(unitarization concept), and 3) the parton radiation effects (hard initial state parton radiation 
in PYTHIA, but not in PHOJET). The comparison of these two models suggests that inferences 
on the parton content in the photon from the measured jet cross sections have an uncertainty 
of about a factor of 2 in the low x-y region. 

The dash-dotted line in Fig. ^ indicates a PYTHIA calculation without multiple interactions 
using a different parameterization of the photon structure (LACl ||3^) which has larger gluon 
content at small x-y than the parameterization used for the other calculations in this paper 
(GRV-LO 10]). This calculation also does not describe the measured energy density (Fig. ^) 
and the jet cross section measurement (Fig. ^). The same PYTHIA calculation with the LACl 
set including multiple interactions gives, after re-adjusting of the transverse momentum cut-off 



to pY^^'^ > 2 GeV (compare section 5.1), a jet cross section which is between 1.5 < if^^ < 2 
and above Eij^^ > 7 GeV a factor of two higher than the calculation with the GRV-LO parton 
distribution function in the photon, and factor 2.5 larger than the data. 

The distribution of the energy density outside the jets in the central rapidity region of the 
jp collision —1 < ??* < 1, i.e. 1 < ij < 3, was used to correct the jet cross section as a function 
of pseudo-rapidity above r/-'^* > 1 for effects of the underlying event: the difference between the 
transverse energy density in the data and the PYTHIA calculation without multiple interactions 
(Fig. 1^: dotted histogram) indicates the event energy outside of jets which is produced beyond 
the hard parton scattering process, its parton showers and their fragmentation effects. Using 
the information on the true parton momentum in generated events the pseudo-rapidity of 
the jet is correlated with xi^^^'^. With the assumption that the calculations of the PYTHIA 
and HERWIG generators represent a good approximation to the energy density outside the 
jets arising from the hard parton scattering process, the transverse jet energies are lowered by 
the energy difference between the data and the calculation without multiple interactions. The 
corrections lower the cross section for jets between 1 < rf^^ < 2 above > 7 GeV by 40%, 
> 11 GeV by 30%, E^^ > 15 GeV by 15%, and demonstrate the strong influence of the 
underlying event energy on the jet cross sections (see also table 0). 



Summary: Properties of jets and jet cross sections 

The width of jets in photoproduction events is observed to decrease with increasing jet energy. 
Comparisons of jet width, measured in 'jp and pp scattering, show that these jets are of a 
common nature. 

Jet cross sections are measured as a function of the jet transverse energy and pseudo-rapidity. 
The shapes of the latter measurements are better described by LO QCD calculations which in- 
clude, in addition to the matrix elements of the parton scattering processes and parton showers, 
interactions of the beam remnants, than by calculations without such additional interactions. 
The strong influence of the underlying event energy on the measured jet cross sections is demon- 
strated. This gives rise to uncertainties in conclusions concerning the parton scattering processes. 

Table 1: Measured differential ep cross-section da/dEx integrated over —1.0 < rj < 2.5 in the kine- 
matical range defined by < 0.01 GeV^ and 0.3 < y < 0.7. 



Et [GeV] 


da/dEr [nh/GeV] 


stat. error 


syst. error 


20-22 


8.1 


0.30 


1.7 


22-24 


5.7 


0.27 


1.3 


24-26 


3.6 


0.19 


0.94 


26-28 


2.3 


0.17 


0.65 


28-30 


1.4 


0.12 


0.44 


30-32 


1.0 


0.10 


0.33 


32-34 


0.80 


0.10 


0.28 


34-36 


0.55 


0.096 


0.21 


36-38 


0.37 


0.065 


0.16 


38-41 


0.26 


0.038 


0.11 


41-44 


0.068 


0.016 


0.032 


44-50 


0.053 


0.020 


0.028 



Table 2: Measured differential ep cross-section da/dE^^ for inclusive jet production integrated over 
-1.0 < 7]^""* < 2.0 in the kinematical range defined by < 0.01 GeV^ and 0.25 <y < 0.7. 



E^' \GeV] 


da/dEf [nh/GcV] 


fildl .error 


sjjsl .vii-or- 


sysl .ovrr-all error 


7.0- 9.0 


20.8 


0.56 


4.4 


5.4 


9.0- 11.0 


5.5 


0.26 


1.2 


1.4 


11.0- 13.0 


1.74 


0.12 


0.37 


0.45 


13.0- 15.0 


0.803 


0.079 


0.17 


0.21 


15.0- 17.0 


0.283 


0.044 


0.059 


0.074 


17.0- 19.0 


0.153 


0.036 


0.032 


0.040 


19.0- 21.0 


0.0715 


0.021 


0.015 


0.019 


21.0- 25.0 


0.0223 


0.0081 


0.0047 


0.0058 


25.0- 29.0 


0.0172 


0.0063 


0.0036 


0.0045 



Table 3: Measured differential ep cross-section da/dE^* for inclusive jet production integrated over 
-1.0 < rf"^ < 1.0 in the kinematical range defined by < 0.01 GeV'^ and 0.25 <y < 0.7. 



E^' \GeV] 


(la/dEf [nh/GrV] 


sinl .error 


■'ijjsl .error 


■'iij.'il .ore rail error 


7.0- 9.0 


10.3 


0.39 


2.2 


2.7 


9.0- 11.0 


2.76 


0.19 


0.58 


0.72 


11.0- 13.0 


0.79 


0.08 


0.17 


0.20 


13.0- 15.0 


0.322 


0.049 


0.067 


0.084 


15.0- 17.0 


0.141 


0.030 


0.030 


0.037 


17.0- 19.0 


0.0672 


0.021 


0.014 


0.017 


19.0- 21.0 


0.0263 


0.012 


0.0055 


0.0068 



Table 4: Measured differential ep cross-section da/dr]^'^* for inclusive jet production integrated over 
E^* > 7.0 GeV in the kinematical range defined by Q'^ < 0.01 GeV'^ and 0.25 < y < 0.7. 





da/dri^^' [nb] 


stat.error 


syst.error 


syst.overall error 


-1.0- -0.5 


5.9 


0.64 


1.3 


1.5 


-0.5- 0.0 


13.2 


0.86 


2.8 


3.4 


0.0- 0.5 


17.8 


1.1 


3.7 


4.6 


0.5 - 1.0 


21.6 


1.1 


4.5 


5.6 


1.0- 1.5 


24.8 


1.2 


5.2 


6.5 


1.5- 2.0 


30.1 


1.3 


6.3 


7.8 



Table 5: Measured differential ep cross-section da/drj^'^* for inclusive jet production integrated over 
Ei^* > 11.0 GeV in the kinematical range defined by Q"^ < 0.01 GeV"^ and 0.25 <y < 0.7. 





da/drf^^ [nh] 


stat. error 


sy.st. error 


syst. overall error 


-0.5 - 0.0 


1.23 


0.20 


0.26 


0.32 


0.0 - 0.5 


1.85 


0.24 


0.39 


0.48 


0.5 - 1.0 


2.36 


0.30 


0.50 


0.61 


1.0- 1.5 


3.30 


0.35 


0.69 


0.86 


1.5- 2.0 


3.27 


0.33 


0.69 


0.85 



Table 6: Measured differential ep cross-section da/drj^'^* for inclusive jet production integrated over 
E:^ > 15.0 GeV in the kinematical range defined by < 0.01 GeV^ and 0.25 <y < 0.7. 





da/dri^^' [nb] 


stat. error 


syst. error 


syst. overall error 


0.0 - 0.5 


0.50 


0.12 


0.11 


0.13 


0.5- 1.0 


0.48 


0.12 


0.10 


0.13 


1.0 1.5 


0.G7 


0.15 


0.14 


0.17 


1.5 - 2.0 


0.65 


0.14 


0.14 


0.17 



Table 7: The measured inclusive differential jet ep cross-section da/drj^'^* after a subtraction of that 
part of the underlying event energy which is found above the calculated energy of the PYTHIA QCD 
generator without interactions of the beam remnants. The error bars reflect the statistical and systematic 
errors, added in quadrature. The kinematical range is defined by < 0.01 GeV"^ and 0.25 < y < 0.7. 



El"^ threshold 


du/drf"^ [nh] 


du/drf"^ [nh] 


[GeV] 


for 1 < rf""^ < 1.5 


for 1.5 < r]^"^ < 2 


7 


14.1 ±6.9 


17.8 ±8.7 


11 


2.2 ± 1.0 


2.2 ±0.9 


15 


0.54 ± 0.25 


0.52 ± 0.24 



7 Conclusions 



In this paper detailed studies of high Et photoproduction events with center-of-mass energies 
around ^s^p = 200 GeV were presented. The selected events have large transverse energy 
Et > 20 GeV in the hadronic final state or exhibit the production of jets with transverse 
energy above E^* > 7 GeV summed in a cone of size R = 1. 

1. ) The differential cross section as a function of transverse energy was measured between 

20 < Et < 50 GeV in the pseudo-rapidity range -2.5 < rj* < 1 (Fig. |2|a, Tab. I). 
The shape is equally well described by a power law E^" with n = 5.9 it 0.1, and by an 
exponential function e~^^'^ with A = 0.21 it 0.01. 

2. ) The transverse energy flow as a function of pseudo-rapidity is distinct from being flat and 

increases towards the origin of the •jp center-of-mass system (Fig. §b). 

3. ) Both the multi-jet rate and the transverse energy contained in jets increase with the total 

transverse energy (Fig. |3|). 

QCD calculations which include, in addition to the hard parton scattering, interactions between 
the beam remnants, provide better descriptions of these transverse energy distributions than 
models without such additional interactions. 

4. ) The transverse energy density in the central region of the 7p collision outside jets in resolved 

photon interactions reaches 1 GeV/rad per unit area in pseudo-rapidity and ip space and 
is far higher than the energy density observed in direct photon processes, 0.3 GeV/rad 
(Fig. I). 

5. ) The transverse energy depositions show short range correlations and long range anti- 

correlations relative to the transverse energy measured at the origin of the center-of- 
mass system (Fig. ^). 

QCD models without interactions of the beam remnants provide too small underlying event 
energy and show too strong energy correlations. In QCD models with beam remnant interactions 
the additional energy from the secondary interaction is essentially uncorrelated with the primary 
hard parton interaction. Such models not only provide the large underlying transverse energy 
density observed in the data, but also reproduce the energy correlations measured in the data. 
This supports the models including interactions between the beam remnants, although it is not 
inconceivable that a more complete QCD calculation of the primary scattering process may 
reproduce the same effects. 

6. ) The jets observed in ^p collisions become denser around the jet axis with increasing jet 

energy (Fig. [pa). 

7. ) The jet cross section as a function of jet transverse energy decreases as (i?jf*)~" with 

n = 6.1 lb 0.5 between 7 < E^* < 29 GeV in the pseudo-rapidity range — 1 < r/-'^* < 2 
(Fig. ga, Tab. |, |). 

8. ) The jet cross section was measured as a function of pseudo-rapidity (Fig. 0b, Tab. ^, P). 

The strong influence of the treatment of the underlying event energy on the cross section 
was demonstrated (Tab. |^). 



In order to draw conclusions on the parton scattering process, QCD jet calculations have to be 
compared not only to the jet cross sections (items 7,8), but also at least to the energy density 
of the underlying event (item 4) . 

Comparisons of the 7p results 1,2,4,6,7 with distributions from pp collisions show that not 
only are the jets of similar shape, but also the features of resolved 7p interactions are similar 
to those of pp scattering. Furthermore, the analysis of direct and resolved photon interactions 
enabled, for the first time, the study of the underlying event energy depending on the energy 
available to the photon remnant, revealing new information on the multi-parton dynamics. 

Acknowledgments 

We are indebted to J.Chyla, R.Engel and T.Sjostrand for help concerning the Monte Carlo event 
generators used in this analysis. We are grateful to the HERA machine group whose outstanding 
efforts made this experiment possible. We appreciate the immense effort of the engineers and 
technicians who constructed and maintain the detector. We thank the funding agencies for their 
financial support of the experiment. We wish to thank the DESY directorate for the hospitality 
extended to the non-DESY members of the collaboration. 

References 

[1] W.J. Stirhng, "Proceedings of the HERA Workshop", ed. R.D. Peccei, Hamburg 1987 and 

references therein 

[2] HI Collaboration, I. Abt et al., Phys. Lett. B314 (1993) 436 
[3] ZEUS Collaboration, M. Derrick et al., Phys. Lett. B342 (1995) 417 
[4] HI Collaboration, T. Ahmed et al., Nucl. Phys. B445 (1995) 195 
[5] G. Knies, J. Phys. G, Nucl. Part. Phys. 19 (1993) 1523 

[6] S.D. Ellis, Proceedings "QCD and High Energy Hadronic Interactions", 28th Rencontre de 
Moriond, Ed. J. Tran Thanh Van, Editions Frontieres (1993) 235 

[7] T. Sjostrand and M. van Zijl, Phys. Rev. D36 (1987) 2019 

[8] G. Schuler and T. Sjostrand, Nucl. Phys. B407 (1993) 539 

[9] HI Collaboration, I. Abt et al., DESY 93-103 (1993) 
[10] HI Calorimeter Group, B. Andrieu et al., Nucl. Instr. and Meth. A336 (1993) 460 
[11] HI Calorimeter Group, B. Andrieu et al., Nucl. Instr. and Meth. A350 (1994) 57 
[12] HI Calorimeter Group, B. Andrieu et al., Nucl. Instr. and Meth. A336 (1993) 499 
[13] J.E. Huth et al., FERMILAB-CONF-90/249-E (1990) 
[14] M. Gliick, E. Reya and A. Vogt, Z. Phys. C53 (1992) 127 
[15] M. Gliick, E. Reya and A. Vogt, Phys. Rev. D46 (1992) 1973 
[16] T. Sjostrand, CERN-TH-6488 (1992), Comput. Phys. Commun. 82 (1994) 74 
[17] T. Sjostrand and M. Bengtsson, Comput. Phys. Commun. 43 (1987) 367 



[18] AFS Collaboration, T. Akesson et al., Z. Phys. C34 (1987) 163 

[19] UA2 Collaboration, J. Alitti et al., Phys. Lett. B268 (1991) 145 

[20] CDF Collaboration, F. Abe et al., Phys. Rev. D47 (1993) 4857 

[21] R. Engel, Z. Phys. C66 (1995) 203 

[22] A. Capella et al., Phys. Rep. 236 (1994) 227 

[23] P. Aurenche et al., Phys. Rev. D45 (1992) 92 

F.W. Bopp et al., Comp. Phys. Commun. 83 (1994) 107 

[24] A. Capella et al., Phys. Rev. Lett. 58 (1987) 2015 

[25] R. Engel, F.W. Bopp, D. Pertermann and J. Ranft, Phys. Rev. D46 (1992) 5192 

F. W. Bopp, R. Engel, D. Pertermann and J. Ranft, Phys. Rev. D49 (1994) 3236 

[26] G. Marchcsini ct al., Comput. Phys. Commun. 67 (1992) 465 

[27] G. Marchcsini and B.R. Webber, Nucl. Phys. B310 (1988) 461 

G. Marchcsini and B.R. Webber, Nucl. Phys. B238 (1984) 1 

[28] B.R. Webber, Nucl. Phys. B238 (1984) 492 

G. Marchcsini and B.R. Webber, Phys. Rev. D38 (1988) 3419 

[29] J.M. Butterworth et al., CERN-TH-95-83 (1995) 

[30] ZEUS Collaboration, M. Derrick et al., Phys. Lett. B297 (1992) 404 

[31] AFS Collaboration, T. Akesson et al., Phys. Lett. B128 (1983) 354 

[32] P. Capiluppi et al., Nucl. Phys. B70 (1974) 1 

BS Collaboration, B. Alper et al., Nucl. Phys. BlOO (1975) 237 

[33] UAl Collaboration, C. Albajar ct al., Z. Phys. C36 (1987) 33 

[34] CDF Collaboration, F. Abe et al., Phys. Rev. D45 (1992) 2249 

[35] S.D. Ellis, J. Phys. G, Nucl. Part. Phys. 17 (1991) 1552 

[36] HI Collaboration, S. Aid et al., DESY 95-156 (1995) 

[37] H. Abramowicz, K. Charchula and A. Levy, Phys. Lett. B269 (1991) 458 

[38] G. Kramer and B.C. Salcsch, Phys. Lett. B317 (1993) 218 
S.G. Salesch, DESY 93-196 (1993) 

G. Kramer and S.G. Salesch, Phys. Lett. B333 (1994) 519 
[39] UAl Collaboration, C. Albajar et al., Nucl. Phys. B309 (1988) 405 
[40] UAl Collaboration, G. Arnison et al., Phys. Lett. B132 (1983) 214 



